Continued fraction approximation for the nuclear matter response function 
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We use a continued fraction approximation to calculate the RPA response function of nuclear 
matter. The convergence of the approximation is assessed by comparing with the numerically exact 
response function obtained with a typical effective finite-range interaction used in nuclear physics. 
It is shown that just the first order term of the expansion can give reliable results at densities up to 
the saturation density value. 
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There are many physical issues that require the knowledge of the response function of a medium to an external probe. 
Well-known examples are the electron scattering by nuclei or the propagation of neutrinos in nuclear matter. In order 
to develop a microscopic theory of response functions in finite nuclear systems one usually starts by considering the 
limiting case of an infinite medium. Infinite nuclear matter as a homogeneous medium made of interacting nucleons 
is a very useful and broadly used concept because of its relative simplicity and its connection with the bulk part of 
atomic nuclei. A popular approach consists in using an effective nucleon-nucleon interaction adjusted to describe the 
nuclear matter properties in a mean field approximation. Then, this microscopic description can be extended to finite 
nuclei. 

In a mean field framework the nuclear response functions must take into account the effects of long-range correlations 
by the Random Phase Approximation (RPA) which is the small amplitude limit of a time-dependent mean field 
approach. This is well suited for those excitations which correspond to small amplitude vibrations, the most typical of 
which being the giant resonances and the low-lying collective states For the theory to be consistent, it is necessary 
that the same effective nucleon-nucleon interaction generates the self-consistent Hartree-Fock (HF) mean field and 
C^- \ the RPA correlations which lead to the excitations of the system. 

There are two types of interactions widely used in non-relativistic approaches, the zero-range Skyrme-type forces [2] 
and the finite-range Gogny-type forces [3j]. Skyrme forces are very often used because of their relatively simple 
analytic form which allows for quite complete RPA calculations in nuclear matter Q as well as in finite nuclei 
On the other hand, finite-range forces require heavier computational efforts to calculate RPA responses in nuclei ^_ 
Furthermore, the only existing methods in this case consist in diagonalizing large size matrices in configuration space. 
It would be useful to have alternative methods such as a direct calculation in coordinate space or momentum space of 
RPA response functions, to avoid the increasingly large configuration spaces. This is possible with Skyrme forces [§] 
but in the case of finite range forces the exchange interactions complicate the problem. 

In this work we study an approximation based on a continued fraction expansion of the response function. Our 
aim is to explore a calculational scheme which can be checked in infinite matter and which offers prospects for RPA 
calculations with finite range forces in nuclei. The continued fraction method is known in the literature [hJ and it 
has been used by many authors to study response functions in the quasi-elastic regime (see Ref. [n| and references 
therein). However, it is difficult to know where to truncate the continued fraction expansion to obtain a desired 
accuracy. It is possible to calculate response functions in infinite matter by performing multipole expansions of the 
interaction and to have numerically accurate results [13] to evaluate various approximation schemes. Therefore, the 
present study aims at assessing the speed of convergence of the continued fraction expansion applied to the response 
functions in nuclear matter, using as an example a Gogny force D 1 0] . We show that this expansion gives good results 
as compared with the numerically exact calculations, even at lowest order. 

In Sec. II we recall the basic features of the continued fraction method applied to the determination of RPA response 
functions in an infinite medium, and we show analytically that it gives the correct result in the special case of a 
Landau-Migdal interaction. In Sec. Ill we discuss the results obtained with a finite range interaction of Gogny type. 
Conclusions are drawn in Sec. IV . 
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II. FORMALISM 



A. General framework 



A general two-body interaction in momentum representation depends at most on 4 momenta. Because of momentum 
conservation there are actually 3 independent momenta, in the case of a translationally invariant interaction. For the 
particle-hole (p-h) case we choose these independent variables to be the initial (final) momentum ki (k2) of the hole 
and the external momentum transfer q. We follow the notations of Ref. [l2j and we denote by a = (S, M; T, Q) the 
spin and isospin p-h channels with 5=0 (1) for the non spin-flip (spin-flip) channel, T=0 (1) the isoscalar (isovector) 
channel, M and Q being the third components of 5 and T. The matrix element of the general p-h interaction including 
exchange can be written as: 

t£' a ' ) (q,k lj k 2 ) = 

(q + k 1 ,kr 1 ,(a)|^|q + k 2 ,k 2 - 1 ,( a ')) • (1) 

To calculate the response of a homogeneous medium to an external field it is convenient to introduce the Green's 
function, or retarded p-h propagator G^(q,cJ,ki). From now on we choose the z axis along the direction of q. In 
the HF approximation the p-h Green's function is the free retarded p-h propagator [I3 |: 

r (auk) ibM m 

CzHF(q,^,kiJ - — ■ — -— yj- ■ — , (2) 

lu + e(fci) - eflki + q|) + 177 

where e(k) is the HF single-particle energy corresponding to momentum k, and the Fermi-Dirac distribution / is 
defined for a given temperature T and chemical potential \x as f{k) = [1 + e( c ( fc )~^ t )/ T ]~ 1 . The HF Green's function 
Ghf does not depend on the spin-isospin channel a. 

To go beyond the HF mean field approximation one takes into account the long-range type of correlations by 
re-summing a class of p-h diagrams. One thus obtains the well-known RPA [l3| whose correlated Green's function 
G^p A (q, u), ki) satisfies the Bethe-Salpeter equation: 



GRp A (q^ k i) = G H F(q,w,ki) +G H F(q,^,ki) j ^%V^' a ) (q,ki,k 2 )G^(q,w,k 2 ) . (3) 

Finally, the response function x (Q> u ) m the infinite medium is related to the p-h Green's function by: 

XRP A (q> w ) =9 J |^G^ A (q,w,ki) , (4) 

where the spin-isospin degeneracy factor g is 4 in symmetric nuclear matter and 2 in pure neutron matter. In the case 
of a system of particles without residual interactions the free response is obtained by calculating Eq. with the HF 
p-h propagator Ghf, thus obtaining the well-known Lindhard function xhf- 



B. Continued fraction approximation 

A direct numerical solution of Eq. (J3j) with a general p-h interaction is possible, as it has been shown in Ref. [12] for 
the Gogny interaction. However, such a method is specifically designed for infinite systems and it would be interesting 
to have an alternative method which can be accurate and at the same time can be used in calculations of finite systems. 
We examine now an approximate way to calculate the RPA response function, expressing it as a continued fraction. 
To simplify the writing of the equations we shall employ the following conventions. First of all we omit the variables 
such as q, it) or k as well as indices (a), unless necessary. For instance, equation ((3j) is written as 

Grpa = Ghf + GhfV^iiGrpa • (5) 

Secondly, for any function -F(ki) depending on a momentum ki we denote by (F) its integrated value over momentum 
space. For example, 



(Ghf) = J 7^p- G H F(ki) , 



(6) 
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so that Eq. Q is simply written as 

Xrpa = 3(G RP a) • (7) 
The Bethe-Salpeter equation is an integral equation which can, in principle, be solved iteratively 

Grpa = Grf + GnF^phGnF + GHFVphGHpVphGHF + ■ ■ ■ (8) 
Correspondingly, the RPA response function is written as 

Xrpa = Xhf + .g(G HF (l)VWl, 2)G HF (2)) + s<G H F(l)Vph(l, 2)G H F(2)^ P h(2, 3)G HF (3)) + • • • (9) 

The brackets imply integrations over chains of variables as shown here. 

In Ref. an approximation was suggested by defining an effective interaction V e g(<l, u, T) such that the RPA 
response function is written as 

Xhf , 1fV . 
Xrpa = T- r ■ (10) 

1 - KjffXHF 

In the RPA neglecting exchange (the ring approximation) the effective interaction does not depend on the hole 
momenta ki and k2 so that Eq. l[T0j) is exact if one replaces V c s by Vph. However, it is important to treat direct and 
exchange terms on equal footing, since they are in general of the same order of magnitude. Here, our point of view 
differs from other works where the direct and exchange interactions are treated on different approximation levels [Til ] . 
We express the effective interaction as a continued fraction 

v oS = ^ . (11) 

l _ ^2 Xhf 

" _ VjiXHF 
1 - . . . 

Each term Vi entering this definition is deduced by expanding formally Eqs. fTTjl and (fTUj) in powers of ViXHF and 
identifying with Eq. (J9j) . The explicit expression for the first two terms are: 

5(Ghf^p1iGhf) 

V2 = — vTj^r Vl - (12) 

First, one can notice that the quantities Grf and xhf are complex functions of q, u> and T, and so are the Vi 
and the effective interaction V c g. Second, the calculations of the V, in the infinite medium involve only products of 
functions, which is somewhat easier numerically than the full calculations of response functions where one needs to 
perform matrix inversions [T3|. Third, one see that V\ is just the average of the full p-h interaction over the squared 
free p-h Green's function. Therefore, the continued fraction approximation could be quite useful for calculating RPA 
response functions if one checks how accurate it can be for a general interaction like the Gogny force. This is what 
we shall examine in Sec. III. 



C. An analytical case: the Landau-Migdal interaction 

The convergence of the approximation can be explicitly seen in the schematic case of a p-h interaction of the 
Landau-Migdal form containing I = and i = 1 terms: 

^ P h = <?{/o + /icos0i 2 } (13) 

where for brevity the same notation /j is used for the Landau parameters in the four spin-isospin channels. For such 
an interaction the RPA response function can be analytically calculated (see e.g. Ref. [j): 

XRPA = -, — 2 — r , (14) 
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where v = u>m*/(qkF), F\ = fiNo is the dimensionless Landau parameter and No = gm* hp / (^tt 2 ) is the level density 
at the Fermi surface, with m* being the effective mass. 

To compare with the continued fraction approximation, we have to evaluate V e g using the interaction ifl3|) . It is 
sufficient to write explicitly the first 3 terms of the expansion of V e s and to obtain the complete series by recursion. 
The integrations involving Ghf have to be carried out in the Landau limit, i.e. q = 0, but finite v. We get: 

Vi = ./o + /i^ 2 (15) 

V 2 = -- — 16 

3 ViXHF 

_ 1 /ii> 2 

V3 ~9WA~ V2 (17) 

It is worth noticing that direct and exchange terms have been treated on the same footing in calculating the V^s. Of 
course for /i = only V\ is needed and one gets the exact result. The effective interaction is 




(18) 
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One can see that this V g leads to the exact result (fl4|l for the RPA response function. 



III. RESULTS FOR A GOGNY INTERACTION 



In this section we apply the continuous fraction method to calculate response functions in infinite symmetric 
matter for a realistic case, using the Gogny effective interaction Dl 0. We choose this parametrization because at 
the mean field level there is a compensation between the direct and the density-dependent contributions. Thus, it 
may be expected that the relative contribution of the exchange term will be somehow enhanced. The purpose is to 
demonstrate the feasibility and rapid convergence of the method. We only present results at T = 0, for which the 
effects of the residual interaction are stronger. 

The task of calculating the V^'s involves carrying out integrals over an increasing number of variables. We find 
convenient to use a multipole expansion of both the HF propagator Ghf and the p-h interaction V p h, as we did in 
the numerically exact calculation of Ref. [T3 |: 

Ghf(<Zj ki) = 5^G/(g,w,fci)y«)(l), 

I 

V ph (qMM) = 5Z«/(g > * 1> fej)y^(l)yi m (2) . (19) 

This allows to get rid of all integrations over angles and we are left with only integrals over the absolute values of 
momenta. For instance, we have 

W = (^I>«G ( >, (20, 

where the integrals implicit in the brackets refer now to the moduli ki. Similar expressions can be obtained for other 

Vi's. 

We can have an idea of the convergence rate by comparing the functions Vixhf and V^xhf- This is shown in Fig.Q] 
for the case of a momentum transfer q=27 MeV. Notice that the scale used to plot V2XHF is about a factor of ten 
larger than that of ViXhf- It can be seen that the imaginary parts of V2XHF are close to zero for the four spin-isospin 
channels. The real parts are generally small compared to 1, but the situation seems less favorable in the channel 
{S,T) = (0,0). From the behavior shown in Fig.l one can expect a rapid convergence of the calculated responses 
already at the level of V2, although perhaps slower in the case of the (0, 0) channel. 

We now examine the strength functions 



S(q,u>) = lmx{q,uj) 



(21) 
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FIG. 1: (Color online) Real (solid line) and imaginary (dashed line) part of ViXhf (top row) and V2XHF (bottom row) for Dl 
interaction in nuclear matter at saturation density po- The transferred momentum is q=27 MeV. The (S, T) channels are shown 
in each panel. 



obtained at various levels of approximation, as compared with the direct numerical solution of Eq. ((3|) presented 
in Ref. [13]. In Figs. [2][3] we show the RPA strength functions for two values of the momentum transfer, at about 
/cf/10 and fcp. The first order gives a reliable description of the strentgh function for all channels except (0,0) as 
expected from the previous analysis. For the (0, 0) channel it is necessary to include the second order. Notice that the 
agreement is independent of the value of q, as no expansion in powers of q has been done. Indeed, as it can be seen 
in Eq. ifTTj) the convergence of the approximation for the effective interaction does not rely on q but on the functions 

ViXHF- 
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PIG. 2: (Color online) RPA Strength function (open diamonds) compared with continuous fraction approximation (1st order: 
dotted line, 2nd order: dashed line) calculated with Gogny Dl interaction in symmetric nuclear matter, at saturation density 
po and momentum transfer q=27 MeV. The thin lines represent the uncorrelated HF strengths. 

However, as the density increases the convergence is deteriorating. In Fig.[4]are plotted the strength functions S^ 0,0 ) 
and the functions ViXRF, V2XHF at density p = 2po in the (0, 0) spin-isospin channel. It can be seen that V2XHF is no 
longer small as compared to VxXhf and consequently, a reliable strength function should require at least the inclusion 
of third order terms in the effective interaction. On the other hand, for densities smaller than po the approximation 
V c ff = Vixbf is sufficient to get accurate results. Of course, the specific convergence found in each channel (S,T) 
depends on the specific interaction used. 

Let us remind that the present approximation is not related to the relative importance of the direct and exchange 
contributions to the particle-hole interaction. Had the exchange term be small as compared to the direct one, a good 
approximation for the response function could be obtained by treating exactly the contribution of the latter and using 
some approximation for the contribution of the former term. This is not the case for the Dl interaction, as it can 
be seen in Table 1. The explicit expressions of these terms are given in Ref. [12] • As the exchange term depends on 
momenta k\, &2 and their relative angle, in the Table are plotted the monopole contributions at the Fermi surface 
(fci = ki = Icf), for two values of q. 
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PIG. 4: (Color online) Same as Figs. [T] and [2] for the channel (0,0) and for p = 2po. 



IV. CONCLUSIONS 



We have examined the efficiency of the continuous fraction method for calculating RPA response functions in infinite 
nuclear matter using a typical finite range effective force. This issue originates from the need of having self-consistent 
theoretical predictions of nuclear responses calculated with realistic interactions. 

We have found that, with the Gogny interaction Dl the continued fraction method is very efficient and the exact 
RPA response functions in the 4 spin-isospin channels are well reproduced already at first order. This is true when the 
nuclear density is of the order of, or less than the saturation density value. At higher densities it becomes necessary to 
include second and higher order terms. The rate of convergence is controlled by the decrease of the terms of successive 
orders VjXhf- In our expansion the direct and exchange interactions are always treated on equal footing. This is 
important since in the nuclear case usually there occurs a strong cancellation of two large numbers, see Table I. 

The encouraging results obtained in infinite nuclear matter open the way to important developments. For example, 
the continuous fraction method for response functions provides a simpler way to evaluate the propagation of neutrinos 
in dense matter such as inside neutron stars. The accuracy of results is under control by the rate of decrease of 
the successive terms ViXHF- In finite nuclei, response functions can be calculated consistently with realistic effective 
interactions without diagonalizing RPA matrices of extremely large dimensions. This can be of some advantage for 
studying heavy and/or deformed nuclei. 
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(S,T) channel 


(0,0) 


(o,D 


(1,0) 


(1,1) 


q (MeV) 


27 


135 


27 


135 


27 


135 


27 


135 




885 


1129 


-363 


-459 


845 


798 


-46 


-146 




-1147 


-1147 


917 


917 


-420 


-420 


583 


583 



TABLE I: Direct (D) and exchange terms (E) in MeV.fm 3 of the Dl p-h interaction in nuclear matter, for p = po and angular 
momentum £=0. 
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